NonOhmic Conduction as a Probe of Coexisting Phases in Manganites 
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A new scaling formalism is used to analyze nonlinear I-V data in the vicinity of metal-insulator transi- 
tions (MIT) in five manganite systems. An exponent, called nonlinearity exponent, and an onset field for 
nonlinearity, both characteristic of the system under study, are obtained from the analysis. The onset field 
is found to have an anomalously low value corroborating the theoretically predicted electronically soft 
phases. The scaling functions above and below the MIT of a polycrystalline sample are found to be same 
but with different exponents which are attributed to the distribution of the MIT temperatures. The exis- 
tence of field scales in manganites underlines the universality of the assumed scaling form in disordered 
systems. 

PACS numbers: 75.47.Lx, 71.30.+h, 72.20.Ht 
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In an intermediate doped manganite A a .Bi_ x Mn03 
(A= trivalent rare -earth and B= divalent alkaline-earth), a 
metal-insulator transition (MIT) usually manifests itself in 
the form of a resistance peak at a temperature Tmi close to 
the Curie temperature T c and is accompanied by colossal 
magnetoresistance. The resistance curve is rather asym- 
metric with a sharp drop on the lower side of T MI , espe- 
cially in single crystals. Spatial microscopy experiments 111 
2] and theories Q 0] attribute this to nanoscale phase sep- 
aration near the transition where the spin-sensitive trans- 
port occurs through a percolative network of ferromag- 
netic metallic (FM) and paramagnetic insulating (PI) do- 
mains. Features such as formation of a backbone of FM 
domainslH ID] at Tmi, significant increase in noise upon 
approaching Tm.t Jal indeed conform qualitatively to the 
standard percolation model H. However, the nature of the 
MIT which is first-order and large deviations in transport 
exponents in manganites from those in standard percola- 
tion model indicate that the present understanding of the 
coexistent phases is still incomplete fl. 

The intrinsic inhomogeneity due to phase separation 
creates conditions favorable for nonlinear (or nonOhmic) 
transport due to wide variations of local electric fields in- 
side a sample. Study of field dependent transport in man- 
ganites through the phase separated regions can be used as 
a tool to probe further complexities in conduction mecha- 
nisms beyond those revealed in zero field measurements. 
Nonlinear conduction has been earlier studied in mangan- 
ites mostly in charge ordered (CO) regimes|9|], in materials 
without CO[ 10], and in orbital ordered samples! 1 !]■ Anal- 
ysis of the nonlinear data so far has been devoid of any 
systematics and mostly empirical. 

Nonlinear studies! 12] in composites of metals and insu- 
lators at room temperature, which are supposed to be pro- 
totypes of the percolation model, reveal that above the per- 
colation threshold (below T MI in the present case) and for 
not too large bias, the (electric) field-dependent conductiv- 
ity a(p, F) of a sample is given by the model-independent 



scaling relation 



a(jp,F) 



(1) 



a(p,0) 

Here $ is a scaling function, p is the metal fraction and F is 
the applied dc electric field. The field scale F (p) is given 
byH 

F (p)~a a *>, (2) 

where <r {p) = a(p, 0) is the zero bias conductivity at 
a metal fraction p. The nonlinearity exponent x p in 3D 
systems was found to have experimental values of about 
0.45 ll2lll3n . The latter is consistent with a model predic- 
tion x p < v jt. Here, v and t are correlation and conduc- 
tivity exponents [7] respectively. In contrast, literature be- 
low the percolation threshold (above T MI ) is rather scant 
in this case. Eq. [Qdescribes a situation where nonlinearity 
is investigated as a function of disorder (i.e. p). However, 
as seen below, nonlinearity can also be investigated as a 
function of temperature, for example, at a fixed or varying 
disorder with a possibility of different exponent x T . Re- 
cently, a similar scaling framework 111 411 has been proposed 
for disordered systems with localized states. In mangan- 
ites, disorder, aside from spin disorder, arises from ran- 
dom potential fluctuations due to A 3+ and B 2+ ion cores, 
Jahn-Teller distortion and local trapping in ferromagnetic 
regions with non-collinear magnetic order[15]. Mangan- 
ites thus provide unique opportunity to study the applica- 
bility of such a scaling formalism. 

In this letter, we present evidence that states around 
metal-insulator transition peaks indeed obey the scaling 
(Eq. Q]) in both single crystal and polycrystalline samples 
of manganites. This proves that origin of nonlinearity is in- 
trinsic, and not due to extraneous factors such as inelastic 
intergrain tunneling. The scaling functions on both sides 
of a MIT in a polycrystalline sample are found to be iden- 
tical albeit with different nonlinearity exponents i T . The 
electric field scale F a is anomalously found to be of the 
order of 1-10 V/cm, about three orders of magnitude less 
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FIG. 1. (Color online) a) Ohmic resistance Rq vs. temperature 
T of two SSCMO samples - one single crystal (with star) and 
another polycrystalline. The symbols, circle and diamond, corre- 
spond to the data of heating and cooling cycles respectively; b) 
Plot of Axt vs. A7Y/2 for four polycrystalline samples (Table 
I). The inset illustrates using the SSCMO data in the panel a how 
AT^ is defined. See text for further explanation. 

than that in granular composites, corroborating the claim of 
electronically soft phases in manganites llloll . In contrast to 
those in granular materials, nonlinearity exponents in man- 
ganites often seem to possess negative signs. 

An important factor in choosing samples for the present 
study (Table I) was the requirement to have resistance 
changes between the baseline and the peak on either side 
of a MIT as large as possible to extract reliably the nonlin- 
earity exponent (see Eq. O. Accordingly, systems with Sm 
were chosen 111 711 . La .275Pr .35Cao.375Mn03_was cho- 
sen for its robust phase separation property [lj, A sin- 
gle crystal of Sm . 55 (Sro.5Cao.5)o.45Mn0 3 (SSCMO* in 
Table I) was grown using the floating zone technique 111 811 . 
Four different polycrystalline samples (ones in Table I 
without star sign) were prepared by the usual solid-state 
reactions. Two SSCMO samples have slightly different 
chemical composition. Of all the samples, LMO is self- 
doped! 19] while rest are the usual hole-doped. The val- 
ues of T MI for different samples analyzed in this work are 
shown in Table I. I-V measurements were done in a Janis 
cryotip with bar-shaped samples placed on sapphire sub- 
strates. Data were taken under constant current condition. 
Maximum current levels were kept low to minimize Joule 
heating in the samples. Both four probe and two probe 
contacts gave the same results. All measurements reported 
here were done at zero magnetic field. 

Figure 1 shows resistance-temperature (R-T) plots for 
two SSCMO samples - one single crystal (with star) and 
another polycrystalline. The single crystal (SC) shows a 
first-order-like sharp MIT transition for T < T MI whereas 
the transition to FM state in the polycrystalline (PC) sam- 
ple is rather gradual. Such behavior, widely known in liter- 
atures, results from a distribution of Tm t 112011 . The hystere- 
sis between cooling and warming cycles is also illustrative 
of first-order phase transition. Steady state I-V measure - 
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FIG. 2. (Color online) (a) Conductance vs. bias in the insulating 
paramagnetic (PI) phase of a SSCMO single crystal at different 
temperatures (T > Tmi) as indicated, (b) Scaling of the data in 
the panel a. Inset shows bias scale, V D vs. linear conductance, 
E for two experimental paths. See text for details. The number 
represents the slope of the linear fit to the data. 



ments at T < T MI are not feasible in the SC sample as the 
sharp transition is metastable with huge relaxation times 
(> 5 hours). The same problem arises in PC samples dur- 
ing cooling cycle but is absent in the heating cycle. This 
makes PC samples a sort of necessity in the present study 
if measurements are to be carried out on both sides of the 
MIT. The S - V (S = l/R) characteristics of SC in the 
warming cycle at different temperatures ranging from 85 to 
105 K (T > Tmi) are plotted on a log-log scale in the left 
panel of Fig. 2(a). The nonlinear response of conductance 
to the application of bias V is clearly seen in the figure. The 
onset bias V (T) is the bias at which conductance starts de- 
viating from its linear value S (T) at temperature T and is 
defined (arbitrarily) such that T,(V a ) ~ 1.1E . The plot of 
the normalized conductance T,(T, V) /S(T, 0) vs. the nor- 
malized bias V I V a using the data of the PI phase of Fig. 2a 
is shown in Fig. 2b which exhibits almost perfect data col- 
lapse, thus validating Eq. Q] The bias scale V Q appears to 
be independent of temperature and is shown in the inset of 
Fig. 2b (circles) as a function of the Ohmic conductance. 
The nonlinearity exponent x^ 1 is thus zero in the PI phase 
of the SSCMO single crystal. The experimental path in this 
case corresponded to cooling the sample down to a temper- 
ature far below T MI and then gradual heating. To check the 
effect of possible hysteresis, the sample in another experi- 
ment was first cooled from the room temperature to 85 K 
slightly above T MI (i.e. still in the PI phase) where heat- 
ing and cooling cycle coincided (Fig. 1), and then heated. 
The exponent x^ 1 in this case turns out to be -0.31 ±0.06 
(squares in the inset of Fig. 2b). Such sensitivity of the ex- 
ponent to the experimental paths perhaps confirm the role 
of hysteresis. 

The fraction of the FM phase for T > Tmi is sup- 
posed to increase as T is decreasedJH |2(i. Such a sit- 
uation is comparable to a non-percolating granular sam- 
ples with different p. We have analyzed such data at a 
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FIG. 3. (Color online) Plots of onset field F vs. Ohmic con- 
ductivity a for different systems as marked. Circles and squares 
correspond to PI and FM phases respectively. Data of LMO and 
Pt-Si0 2 were shifted to right by the factors as mentioned. 



fixed temperature of three samples of Pt-Si0 2 with differ- 
ent metal fractions II2 ill (supplementary materials). Inter- 
estingly, the nonlinearity exponent x p in this case is found 
to be ~ as in SSCMO*. It is seen however from Table 
I that the same exponent x^ 1 in the polycrystalline sam- 
ples can have non-zero values, both positive and negative 
although the magnitudes are not large. The field scale F a 
in the prevalent theories Ell 12211 is given by the condition 
that the field-activated conduction becomes comparable to 
the temperature-activated one i.e., 



eF w = k B T. 



(3) 



Here w, a relevant length scale, may be the mean sepa- 
ration length between two metallic clusters in a granular 
system 12111 or the average hopping distance in a variable 
range hopping svstem l22ll . If w and consequently, F re- 
main independent of p we have x p = 0. This means that in 
a manganite system even though amount of the FM phase 
may vary with T, the average tunneling distance between 
two FM clusters in PI phase remains fixed. This however 
can not explain any non-zero exponent. Let us now con- 
sider another aspect of the field scale, namely its magni- 
tude which is displayed in Fig. 3 for both SSCMO* and 
Pt-Si0 2 . It is seen that the scale in SSCMO* is about three 
orders of magnitude less than that in Pt-Si0 2 . This can be 
easily understood by noting that a carrier moves from one 
metallic cluster to another which are chemically and struc- 
turally same but distinguished only by spin or charge order- 
ing. Hence, the energy barrier that a carrier must overcome 
must be much less than that in a medium of chemically 
dissimilar materials like Pt and Si0 2 . On the other hand, 
taking the separation distance in a manganite as large as 
100 nm[l] a field of 10 V/cm (F Q in SSCMO*, Fig. 3) cor- 
responds to an anomalously low barrier height of about 1 K 
compared to about 300 K estimated from noise studies ]23J. 
Furthermore, such low value of 1 K contradicts the basis of 



10" 



10" 



10 



T<T M| SSCMO (a) 
77K 



o-o — o oo o 



CP 



85 ,^ 



90 

* — *-*—*" 



4* 



T> S 5 0K J» 

o-o-o-coo^ 



125 

A A — A AA A 

113 



-oo-o 



102 



10 10 10 

V(V) 



10 10 10 

V(V) 



2 -0.3 



x™=0.17±0.01 




0.2 4 
10~* 



(c)a 



T<T,„ TO; 



10* 



10" 

~ 4 -^--»---w 4 * 



10" 



10" 



v/v„ 



10" 



10 



FIG. 4. (Color online) The upper panels show plots of conduc- 
tance vs. bias in a SSCMO polycrystalline sample in the ferro- 
magnetic metallic phase (FM) for T < Tmi (a) and PI phase for 
T > Tmi (b) at temperatures as indicated. The lower panel (c) 
shows the data collapse of the same data of the upper panels. The 
scaled curve of the FM side is shifted for clarity. Inset shows two 
log-log plots of Fq vs. (to data corresponding to PI (open cir- 
cle) and FM sides (open square) respectively. The solid lines are 
linear fits to the data with slopes and 1 respectively. 



nonlinearity in Eq. [3] the temperature range of interest be- 
ing 81-95 K. It may be noted that Eq. [3] is not compatible 
with Eq. |2]which is a typical power-law relation in critical 
phenomena. 

The nonlinear I-V curves could also be caused by Joule 
heating. A simple mean field approach shows that V ~ 



or, x 



Joule 



—0.5 [24]. As none of the exponents 
obtained is close to this value it can be concluded that the 
nonOhmic behavior is an intrinsic feature and not caused 
by heating. Let us now turn to measurements on a host of 
polycrystalline samples (Table I). Since results in various 
samples are basically similar, we consider only the results 
for the polycrystalline SSCMO in Fig. 4 for the purpose of 
comparing results from the single crystal of the same sys- 
tem discussed above. S — V plots at various temperatures 
both below and above T MI are shown in Figs. 4a and 4b 
respectively. Note that in both cases conductance increases 
from bottom to top. Curves in panels a and b look qual- 
itatively similar but exhibit a subtle difference as shown 
below. Both sets of data were found to follow the same 
scaling as in Eq. Q]and actually collapse on the same curve 
(Fig. 4c). But the scaled data in FM phase (T < Tmi) has 
been shifted in the figure for the sake of clarity. V thus 
obtained are shown in log-log plots (open squares for FM 
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TABLE I. Sample parameters and nonlinearity exponents. See text for definitions. 



System 


Abbre- 
viation 


Type 


Tmi 
K 


T\/2 

a / 
K 


X FM 


t pi 


Sm .55(Sro.5Cao.5)o.45MnC>3 


SSCMO* 


Single crystal 


82 









Sm .55Sro.3375Cao.ii25Mn0 3 


SSCMO 


Polycrystal 


94 


85 


0.17 ±0.01 


0.09 ±0.01 


Sm .55Sr .4 5 MnC>3 


SSMO 


Polycrystal 


69.5 


59.7 


-0.23 ±0.01 


-0.14 ±0.01 


Lao.275Pro.35Ca . 3 75MnC>3 


LPCMO 


Polycrystal 


113 


102.6 


-0.14 ±0.02 





La .87(Mn 2 3 )o.i3Mn03 


LMO 


Polycrystal 


155 


62 


-0.83 ±0.01 


-0.16 ±0.01 



phase and open circles for PI phase) in the inset of Fig. 4c. 
The solid lines are power-law fits to Eq. |2]wifh the expo- 
nents as indicated. The exponents in two sides of the MIT 
are different, but are positive and non-zero in contrast to 
zero value in the singe crystal. Moreover, V a is generally 
less in polycrystals than in the single crystal (Fig. 3). This 
is contrary to what will be expected if the origin of nonlin- 
earity were due to intergrain effects. It is seen from Table I 
that | x^ 1 | is generally less than | x^ M \. F as a function 
of conductivity for all polycrystalline samples are shown in 
Fig. 3. 

The FM phase at T < T MI possesses a continuous 
metallic path or a backbone through a sample. Conductiv- 
ity of a backbone is too large to make any contribution to 
nonlinearity. The nonlinear conduction above the threshold 
in standard percolation is thought to arise from tunneling 
across thinly separated conducting bonds dangling off the 
backbone lil2ll . in contrast to tunneling of a carrier from one 
metallic cluster to another across insulating barrier below 
the threshold. Thus, the observation of the same scaling 
function above and below T MI in the polycrystalline sam- 
ples was unexpected. Alternatively, as mentioned earlier, 
there exists a distribution of T MI and hence, a fraction of 
PI clusters which gets progressively reduced as tempera- 
ture is decreased 11201 . This is indirectly supported by the 
fact that the sign as well as the magnitude of the exponents 
s™ are generally very different from ~0.45 found in the 
binary composites above the threshold. Conduction across 
such PI clusters may give rise to same scaling functions 
on both sides of T MI . It is then necessary to explain why 
the exponents on both sides of the transition are different. 
We postulate that the difference is related to the width of 
the distribution of T MI , which is taken to be Gaussian lEoll . 
A measure of the width, denoted by ATx/2, is assumed 
in this work to be the difference between T MI and the 
temperature T\/ 2 such that R (Ti/ 2 ) = R (T MI )/2 (see 
inset of Fig. lb). The absolute difference in exponents, 
Axt = | x PI — x^ M | when plotted against ATi / 2 yielded 
a straight line passing nearly through the origin (Fig. lb). 
Of course, no explanation of the individual exponents, par- 
ticularly negative ones, is available at present. 

In summary, we suggested a model-independent scaling 
formalism to analyze the electric field-dependent conduc- 
tion data in manganites. It was shown how quantities such 
as the field scale or nonlinearity exponents could be used 



to gather further information about the underlying physical 
states. There is a need for theoretical modeling to explain 
great variety of exponents found in this work. The same 
scaling analysis could be applied to other regimes such as 
charge ordered states in manganites. 
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